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Abstract 
Rapid urbanization and urban growth, particularly in the developing worlds, is continuing to 
be one of the crucial issues of global change in affecting the physical dimensions of cities. 
This study proposes a technique to extract urban built-up land features from, Multi Spectral 
Scanner System (MSS-1977), Landsat Thematic Mapper (1998, 2003, and 2014) of part of 
Delhi NCT in India as examples. The study selected two indices, Normalized Difference Built- 
up Index (NDBI), Normalized and Normalized Different Vegetation Index (NDVI) to represent 
three major urban land-use classes, built-up land and vegetation, respectively. The 
relationship between land use/land cover (LULC) change and population shift and their 
effects on the spatiotemporal patterns of urban area were quantitatively examined using an 
integrated approach of remote sensing, geographical information systems (GIS). 
Consequently, the seven bands of an original Landsat image were reduced into three 
thematic-oriented bands derived from above indices. The three new bands were then 
combined to compose a new image. As a result, the spectral signatures of the three urban 
land-use classes are more distinguishable in the new composite image than in the original 
seven-band image as the spectral clusters of the classes are well separated. Therefore, the 
technique is effective and reliable. In addition, the advantages of over NDVI and over NDBI 
in the urban study are also discussed in this paper. Furthermore, in combination with the 
detection of LULC change, an analysis of the spatially differential growth rates for developed 
land area and population size revealed an urban & sub-urban gradient pattern of population 
shifting, as evidenced by a sharp increase in developed land area within the middle sub- 
zones at the urban fringe and the exurban sub-zones beyond the outer traffic ring. 
Consequently, changes in LULC and population shifts resulted in significant variation in the 
spatiotemporal patterns of the urban area. 
Keywords: GIS, Landsat Image, NCT, NDBI, NDVI, Remote sensing. 


Introduction 


Population growth is a major problem in unplanned growth of town =  and_ £cities. 
much urban country like China, India, and Africa Urbanization has been a_ universal and 
etc. Rapid Urbanization causes disorganized and important social and economic phenomenon 
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taking place all around the world. Urbanization 
could be the most powerful and _ visible 
anthropogenic force that has brought about 
fundamental changes in land cover and 
the 


Urbanization and urban growth, especially in 


landscape pattern around globe. 
the developing world, is continuing to be one of 


the crucial issues of global change. 


Background of the study 

Urban sprawl is the scattering of new 
development on isolated tracts, separated from 
other areas by vacant land (Ottensmann, 1977). 
Remote sensing has great potentials in studying 
urban environments and urban/suburban 
landscape when high spatial resolution imagery 
is available (Jensen and Cowen, 1999). Remote 
sensing data are more widely used for the 
analysis of pattern and process rather than 
In 2009, over 3.4 


billion people in the world resided in urban 


Causes or consequences. 


areas aS a result of rapid urban growth. This 
figure is estimated to increase to 6.5 billion by 
2050. This urban growth phenomenon has 
raised challenges for many countries in both the 
Urban 
growth is a complicated process involving the 


developing and developed worlds. 


Spatiotemporal changes of all socio-economic 
and physical components at different scales. 


Problem Identification 

Many cities worldwide are growing very 
rapidly, leaving urban and transport planners 
with continuous challenges in planning for 
livable environments. This rapid growth causes 
complex interrelated urban growth = and 
transportation issues. In most of the developing 
countries, urban development is occurring in a 
haphazard manner without concerning the 
ecological issues. Rapid and unplanned urban 
growth is also threatening the sustainability of 


the cities. In this context, it has been a major 
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challenge to the planners and decision makers 
the 
occurring over urban areas. Since, urban growth 


to identify Spatial-temporal changes 
is a continuous process; it can only be identified 


from analyzing long term spatial-temporal data. 


Aim of the Study 

The present study aims to identify the spatial- 
temporal changes occurring over Delhi National 
Capital Territory region using long term satellite 
data. 


Objective of the study 
> To estimate areas under different land 
use and land cover classes using 
supervised image classification 
technique. 
> To identify the — spatial-temporal 
changes through change’ detection 


technique and detect the pace and 
direction of urban growth using long 
term satellite data. 

> I\Identification of urban growth through 
multi-temporal fractional images and 
comparing it with fractional images of 
hyper spectral satellite data. 

> Mapping the socio-economic changes 
of the study area and analyzing the 
causes of such rapid urban growth in 
Delhi. 


Study Area 

Delhi is the capital of India. The state is 
spread over an area of 1473 square kilometer. 
According to the Indian geography the state is 
the the 
subcontinent, amidst the ranges of Himalaya 


located at center of Indian 


and the Ravalli. 


Location 

The study area is located Centre of the Delhi 
in India. The study area is bounded by Latitude 
28°24'16.84"N to 28°52'52.31"N and Longitude 
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77°5'12.02”. E to 77°10'9.83"E. The total 
boundary of the study area is 1473.1677sq.km. 


LOCATION MAP OF STUDY AREA 


HORVAHA 4 NCT OF DELHI 








TE Tare Traore rr arere 


Fig. 1. Location Map of the Study Area. 





Topography 

The entire topography of Delhi is divided into 
a ridge, the Yamuna Flood Plain, the Plain. It is 
interesting to note here that each of these 
regions is marked by distinct type of 
vegetation. The ridge area of the city offers the 
right factors that favor the growth of acacias 
and other cacti. However, during the monsoon, 
herbaceous plants grow in abundance in the 
ridge. As far as the plain region of Delhi is 


concerned, it is characterized by shish am trees. 


Drainage 
Drainage pattern of an area is very 
important in terms of its groundwater 


potentiality. It is the source of surface water 
and is affected by structural, lithological and 
geomorphologic control of an area. This may be 
due to more or less homogeneous lithology and 
structural controls. Yamuna River is the main 
control of drainage system along the NCT of 
Delhi. 

Reverie type of vegetation grows along the 
plain of Yamuna. National Capital Region in 
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general, is a part of well integrated drainage 
system of the Ganga basin. Drainage is an 
important element of physical infrastructure 
and constitutes removal and disposal of surplus 
rain/irrigation water from the land. 


Vegetation 

Vegetation cover in this area has been found 
to be predominantly. According to report, the 
total vegetation cover is 367.425 sq. Km., 
279.3114sq. Km., 404.3061sq. Km. and 
71.5851.Km. in 1977, 1998, 2003 and 2014 
respectively. The Vegetation of Delhi varies with 
its varied topography and comprises small and 
medium sized plants and shrubs. 


Climate 
Delhi 
climatic conditions that are faced by this city. 
Delhi is 
characterized by scorching heat and unbearable 


Weather varies with the different 


Summer weather condition § of 
temperature. Temperature reaches almost 45°C 
in the summer months. On the contrary, winter, 
which lasts from December to January, is 
extremely cold. Temperature falls to almost 5°C 
during the winter months in Delhi. The rainy 
season in Delhi begins in June and continues 
almost till October. Delhi receives most of its 
rain during this period from the Northwesterly 
winds. Most of the precipitation occurs in the 
month of July. 


Data Collection 
Urban and 


identification land use change in Delhi using 


Detection of growth 


Landsat satellite image. Used temperature 


mapping algorithm for detection urban growth 
and Land use & land cover classification have 
done maximum _ likelihood 


been using 


classification technique. 
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Table 1. Details of the satellite data used in this study. 


Pata ee 


Google Earth Image 2014 Google Earth 


Table 2. Characteristics of the Landsat 7 TM, MSS and Landsat 8 sensor and the using data was 
collected from USGS web site. 


| Sensors Spectral bands Spatial resolution Spectral region Time of acquisition 


Land sat-TM 30M VNIR 1998, 2003 
February 





Band 6: 1.04-12.5um 120 M 
Band 7: 2.09-2.35um SWIR 


Land sat-MSS Band 1:0.50-0.60um 
Band2:0.60-0.704m (Red 
Band 3:0.70-0.80um 


Reflected IR 
Band 4:0.80-1.11um 


February 1977 
Land sat 8 Band 1:0.45-0.52um 2014, 
February 
30m 
-O. 5m 


Band 8: 0.52-0.90um Panchromatic 


Satellite image (Land sat-TM 1998, 2003, & Land sat-8 2014) 
Raster map of the study area 


Delineation of the study area 
(Image based) 


Digital Indexing Band layer stack 


Correlation study 
NDVI 


Built-up Land Extraction 


Geo-Referencing Geometric & Radiometric Correction 











Figure 2. Methodological flow chart. 
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7 \GHye; 
a —— 
Source Landsal TM Feb2oo0d, in Delhi Sours: Lancaat Th Fab A014, in Gelhi 


Figure 4. Normalized Difference Vegetation Index in Delhi, 2003 &2014. 
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Legend | 
== Built-up land a ae ; ~—*" Legend 


: 4: ; is TER RF 
Source:Landaat TM.Feb, ee , im Delhi hg | Built-up land 
Source:Landear TMF eb2003, on Delhi 


Figure 5. Normalized Difference Built-up Index in Figure 6. Normalized Difference Built-up Index in 


Delhi, 1998. Delhi, 2003. 


m Build-up Area 


2014 2003 1998 


Legend Year 


a Built-up land 


SouronLandsat # Feb 2004, in Cadhi 


Fig. 8. Normalized Difference Built-up Index in 


. : . . Delhi, 1998, 2003 & 2014 with Column diagram. 
Figure 7. Normalized Difference Built-up Index 


in Delhi, 2014. 
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Different image processing & enhancement 
technique used for identification the time- 
sequential changes in land use parents. 


Satellite data 
The data used for this dissertation consists 

of- 
(1) Landsat7 MMS imagery, from band 1to band 
4, on February, 1977. 
(2) Landsat7 TM imagery, from band1 to band 7 
on February, 1998, 2003. 
(3) Landsat 8, 2014 imagery, from band1 to 
band 8 on February. 

Here, Land sat TM, MSS, and Landsat 8 
images are used to fulfill this work 

Table-2: Characteristics of the Landsat 7 TM, 
MSS and Landsat 8 sensor and the using data 
was collected from USGS web site. For the 
Extraction of Built-up land | used the MidIR, NIR 
band. Spatial resolution of NIR band of Landsat - 
TM & Landsat-MSS is 120m and 60m 
respectively and spectral resolution of visible, 
infrared & middle infrared region of Land sat- 
™T™ & MSS is 30m & 60m. Due its moderate 
resolution spatial coverage, Land use & land 
cover classification has been done for change 
detection analysis. 


Instruments Used 

Garmin GPS 12: GPS was used to obtain the 
geographical co-ordinates of the observed filed 
location during the ground truth study for 
distributed point, & 
use/land cover information. 


collecting Well land 


Software’s used 

Image Processing: ERDAS Imagine 9.0 & ENVI 
5.0 

Geographical Information System: Arc GIS 10.1 
Others: MS-Office-2010(MS-word, MS-Excel) for 
documentation and calculation. 
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Data analysis 

The primary source of the data used for the 
built-up land investigation on satellite image. 
The remotely sensed data used in this test is a 
Landsat 8, TM/ MSS image. The images are 
cloud-free and have excellent quality. 
Normalized Difference Vegetation Index 
(NDVI) 

Rouse et al., (1974) developed the generic 
Normalized Different Vegetation Index (NDVI). 
In the case of this study, NDVI has been used to 
validate one of the objectives. 

NIR — RED 
NIR+ RED 
A Normalized Difference Vegetation Index 


NDVI = 


(NDVI) is an equation that takes into account 
the amount of infrared reflected by plants. Live 
green plants absorb solar radiation, which they 
use aS a source of energy in the process of 
photosynthesis. The reason NDVI is related to 
vegetation is that healthy vegetation reflects 
very well in the near-infrared part of the 
electromagnetic spectrum. Green leaves have a 
reflectance of 20% or less in the 0.5 to 0.7 
micron range (green to red) and about 60% in 
the 0.7 to 1.3 micron range (near-infrared). This 
spectral reflectance is themselves ratios of the 
reflected over the incoming radiation in each 
spectral band individually; hence, they take on 
values between 0.0 and 1.0. and shows the 
Normalized Difference Index (NDVI) drives from 
Land sat satellites image of the year of 2001 
and 2011. In Figure: 2 NDVI values 0.0814 to 
0.4375 represents the green coverage surface 
and Figure: 3 NDVI value 0.1559 to 0.4516 
represents the green coverage surface. 


Index (NDBI) 
These equation can be used for monitoring 

the built-up land Zha et al., (2003) calculated a 

normalized Difference Built-up Index (NDBI). 


Normalized Difference Built-up 
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_  MidIR — NIR 
NDBI = Bui — NDVI?“ = waar | NIR 

Where, this resulted in an output image that 
contained only built-up and barren pixels having 
positive values while all other land cover had a 
value of O or -254. The technique was reported 
to be 92 percent accurate. 

The built-up areas of 2001 and 2011 were 
extracted from NDBI images. To extract the 
built-up area from the NDBI images we have to 
adopt the thresholding technique, where a 
creating range of DN value is extracted from the 
index raster which may represent some typical 
kind of feature of interest. In general, we know 
that the positive values of NDBI are represent 
the built-up feature and negative values is 
represent the other land surface features. 


Conclusion 

The land-use/land-cover map prepared from 
Landsat-TM, MSS & Landsat 8 data _ for 
identifying the time sequential change of LULC 
patterns. Built-up area has rapidly increased 
during 1977 - 2014. It is also noteworthy that 
the rate of urban growth was highest during 
2003 to 2014 and a total increase of high 
density built-up land was found 47.01%. It was 
happened mainly at the expense of agricultural 
land located at the outskirt of the city. 
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